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Abstract 

We discuss QCD in the Hamiltonian frame work. We treat finite density QCD 
in the strong couphng regime. We present a parton-model inspired regularisation 
scheme to treat the spectrum (^-angles) and distribution functions in QEDi^i. We 
suggest a Monte Carlo method to construct low-dimensionasl effective Hamiltonians. 
Finally, we discuss improvement in Hamiltonian QCD. 



1. Introduction 

There are presently the following problems in QCD: How to compute non- 
perturbatively (on the lattice): (a) QCD at non-zero chemical potential - finite 
quark density, (b) Hadron structure functions, at small xb and small Q. (c) S- 
matrix, scattering cross sections, phase shifts, (d) Decay amplitudes. Here, we 
want to discuss recent progress in the Hamiltonian formulation, QCD at non- 
zero chemical potential in the strong coupling regime, distribution function in 
QEDi^i, the Monte Carlo Hamiltonian, and improvement of the continuum 
behavior of the Hamiltonian. 



2. Finite density QCD 

According to the big bang model of cosmology creation of quarks and gluons 
occured at about 10~^^ sec, quark-gluon plasma {QCP) was formed at about 
10~^^ sec, then there was a phase transition and the formation of protons 
and neutrons began at 10~^ sec. Today the universe is in a low temperature 
and low density phase: quarks are confined. Quark-gluon plasma may exist 
in cores of very dense stars like neutron stars. Experimental search for QCP 
is carried out at the Relativistic Heavy Ion Collider (RHIC) at BNL and at 
the Large Hadron Collider (LHC) at CERN. Theorists search for a precise 
determination of phase structure of QCD at finite temperature T and finite 
chemical potential n (finite density). Shuryak [1] has predicted that QCD 
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at high temperature and density displays beyond the QGP pase transition a 
much richer phase structure, due to strong instanton effects. 

Here we suggest to consider finite temperature and finite density QCD in the 
Hamiltonian formulation. Consider the grand canonical partition function at 
temperature T and chemical potential // 



Z = rre-^(^-''^), N = I d^x7/j''{x)j/j{x), P ^ {keTyK (1) 

H is the Hamiltonian and N is the fermion (quark) particle number operator. 
No complex action problem occurs! Let us consider as example the energy 
density for free quarks 

e - --TrHe-^^"-^'^^ - --^^\ a (2) 



Taking the zero-temperature limit (T — > 0) and the chiral limit (m — > 0) 
the energy density density (with the contribution /i — subtracted) becomes 

2 

^sub — The same result is obtained by starting from a lattice Hamiltonian 
with Wilson fermions. 



Strong coupling QCD at non-zero chemical potential. 

Lattticc QCD at /i = confines quarks and spontaneously breaks chiral sym- 
metry. Note: The strong coupling regime << 1 is not compatible with 
continuum limit a — > 0. The goal of the following calculation is to get a better 
understanding of mechanism of chiral phase transition. We make an ansatz for 
the lattice Hamiltonian in the strong couphng regime following Ref.[2]. Then 
the following results are obtained [3] . The vacuum energy becomes 



En = 2NfN,N, Um^l - /.)e(/. - m^O - - 



(0) 



,(0) 



(0) 
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Cj,^{Nl-l)/{2N,), 



(3) 



where rrS^y^^ is the dynamical quark mass at = 0. In order to compute the 
chiral condensate and the critical chemical potential one can use the Feynman- 
Hellmann theorem. The chiral condensate at = becomes 



<ipip > 



(0)- 



-2Nc[l-4d/{g'Cj,)], 



(4) 



and the chiral condensate at /i ^ becomes 



(5) 
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As a result one has 



< iptl) for // < m^y^ : chiral symmetry broken, 

< >= for /X > mfy^ : chiral symmetry restored. (6) 

Thus we make the prediction of a first order chiral phase transition. The 
critical chemical potential is fXcrit = t^T'dyn- Our result for Hcru is in agreement 
with the result obtained by Le Yaouanac et al.[4] However, Le Yaouanac et 
al. find a phase transition of 2nd order. There have been predictions of the 
nucleon mass [4] M^^^^ ^ 3m^°^„. Thus our result is fjLcrit ~ ^nud/'^i being 
consistent with other theoretical predictions [4] . We have moreover computed 
the quark number density, the corresponding susceptibility and some thermal 
masses [3]. 

3. Hamiltonian in a fast moving frame 

At hand of the massive Schwinger model, Q-E-Di+i, we want to discuss the 
following topics: (i) A physically motivated many-body method to solve a 
Hamiltonian many-body problem and compute structure functions, (ii) The 
treatment of the ^^-angle in Hamiltonian gauge theory. The motivation for the 
Hamiltonian many-body method comes from the parton model of deep inelas- 
tic lepton-hadron scattering. Experimentally, the Breit frame is a convenient 
choice for the interpretation of structure functions. We suggest to use a fast 
moving frame also in a computational study [5]. The Hamiltonian in axial 
gauge = is given by 

I dx^ {Tjj^Hdsij + mTjjTjj) + y / f^^' (V'V) ^ (V'V)- (7) 

-L -L ^ 

Wc start from a space-time lattice of spacing a and length [—L,L]. Then we 
go over to a momentum lattice with cut-off A = n/a and resolution Ap = 
tt/L. We use the assumption (parton model): Left moving particles are not 
dynamically important, when p/i?/sica/ particles are considered from a reference 
frame with velocity v = P/E ^ c. On the momentum lattice the following 
bounds hold: < p^,p^ < P. A fast moving particle is Lorentz contracted. It 
fits into a small lattice volume (compared to rest frame). A small lattice size 
implies a small number of virtual particle pairs created from vacuum: Number 
of vacuum pairs oc vacuum density x lattice size. For the scaling window holds: 
a < Compton wavelength < L and a < < L, where ^ is the correlation 
length in dimensionless units and M = ^ is the physical mass of particle 
in the ground state. In the strong rclativistic regime holds M << P < A. 
In the limit a — > the scaling window is replaced by: 1/P < 1/M < L. 
In Ref.[5] we have computed the mass spectrum: Vector boson and scalar 
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boson, and compared it with predictions of chiral perturbation theory. We 
also determined the dependence on the 9 angle, defined by a global gauge 
transformation ^)[A] e"'" ^'[A], n = 0, ±1, ±2, . . .. The 6'-action S and 
topological charge q are related by 6*] = S[A] — 9q[A]. In Ref. [5] we have 
also studied the massive Schwinger model and its dependence on the ^-angle 
in the chiral limit. For the massless Schwinger model the chiral condensate is 
analytically known, 

= - — j= cos{e) « 0.1599 cos(^). (8) 



Using the Feynman-Hellmann theorem to relate the vector mass and the chiral 
condensate, we determine the slope from numerical data and obtain Ri 
0.16 cos(^), which is close to the analytical result. 

The Hamiltonian approach also allows to compute wave functions, giving in- 
formation on the parton structure. 

fixs) =< ^v{P)\l[blh + 44|*y(P) >, (9) 



where Xb = ^ is the fraction of momentum of vector boson carried by parton. 
One observes that a large lattice (N=200) is needed to resolve distribution 
functions precisely (see Ref. [5]). 



4. Stochastic Hamiltonian 

Lagrangian lattice field theory and its non-perturbative numerical solution by 
Monte Carlo methods is very successful. The typical number of configurations 
(depending on the statistical error) is in the order of ~ 10^. In Ref. [6,7], we 
have propopsed to construct an effective Hamiltonian, by treating the many- 
body problem by applying as far as possible Monte Carlo with importance 
sampling (e.g. Metropolis algorithm) in analogy to Lagrangian lattice simu- 
lations. Let us consider as example 1-D Q.M. We compute transition matrix 
elements in imaginary time 

Mij{T) =< Xile-^'^/^lxj >, i,jel,...,N, (10) 



where Xi are discrete position eigen states. We define an effective Hamiltonian 
via eigenfunctions and eigenvalues: 

M{T)^U'^D{T)U, Ul^<e,\Elff >, D,{T) ^ e'^i'^^'K (11) 



We compute matrix elements using Monte Carlo with importance sampling 
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Fig. 1. Harmonic oscillator. Specific heat C{/3) from Monte Carlo Hamiltonian, using 
stochastic basis. 

Fig. 2. QCDi+i. Fermionic improvement, x = -{i>ip)sub/i9lat^c) versus l/gf^t for 
Nc = 3 with Wilson fermions. Crosses: r = 0.1, Diamonds: r = 1. 

Fig. 3. QCDi^i. Fermionic improvement, x = -{'4^'^)sub/{giatNc) versus l/gfat for 
Nc = 3 with improved Wilson fermions. Crosses: r = 0.1, Diamonds: r = 1. 

(e.g. Metropolis algorithm). Splitting the action S = Sq + Sy, we treat O — 
exp[—Sv/h] as observable and compute the matrix element by Monte Carlo 

Mij{T) = M^^\t) X J[dx] exp[-Sv[x]/h] exp[-So[x]/h]\::fo 

/ j[dx\exp[-SMmi'-l- (12) 

The construction of a stochastic basis has been proposed in Ref.[7]. The result 
of a numerical study is shown in Fig.[l]. 

5. Improvement in Hamiltonian QCD 

Szymanzik [8] suggested improvement program for lattice field theory: Con- 
struction of "improved" action, such that error due to finite lattice spacing a 
becomes smaller. Progress in improvement for lattice pure gauge theory was 
introduced by Lepage et al. [9], suggesting tadpole improvement. We have 
suggested how to construct an improved Hamiltonian on the lattice. For pure 
gauge theory the construction of an improved Hamiltonian for SU{'i) gauge 
theory via the transfer matrix has been reported by by Luo et al. in Ref.[10]. 
Improvement of the fermionic sector of a lattice Hamiltonian has been carried 
out for QCDi^i by Luo et al. in Ref. [11]. A numerical study was carried out by 
Jiang et al. [12]. An example is shown in Fig. [2,3], showing that improvement 
reduces the dependence on Wilson's r-parameter. 

6. Conclusion 

Hamiltonian formulation of QCD is attractive, for the following reasons: (a) 
There is no sign or complex action problem! (b) It is suitable for the com- 
putation of wave functions and structure functions, (c) It is suitable for the 
computation of the S-matrix. The problem is: How to solve the Hamilto- 
nian many-body system? (d) Proposal: Construct an effective Hamiltonian 
via Monte Carlo, (e) Improvement program: Also possible in Hamiltonian for- 
mulation. 
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